Ultrafast electron dynamics in metals 
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During the last decade, significant progress has been achieved in the rapidly growing field of the 
dynamics of hot carriers in metals. Here we present an overview of the recent achievements in the 
theoretical understanding of electron dynamics in metals, and focus on the theoretical description of 
the inelastic lifetime of excited hot electrons. We outline theoretical formulations of the hot-electron 
lifetime that is originated in the inelastic scattering of the excited quasiparticle with occupied states 
below the Fermi level of the solid. First-principles many-body calculations are reviewed. Related 
work and future directions are also addressed. 
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INTRODUCTION 

Recent advances in femtosecond laser technology have 
made possible the investigation of electron transfer pro- 
cesses at solid surfaces, which are known to be the basis 
for many fundamental steps in surface photochemistry 
and ultrafast chemical reactions. Q, |^ |^ These are typi- 
cally atomic and molecular adsorption processes and cat- 
alytic reactions between different chemical species, which 
transfer energy from the reaction complex into the nu- 
clear and electronic degrees of freedom of the solid sub- 
strate. These reactions may induce elementary excita- 
tions, such as quantized lattice vibrations (phonons), col- 
lective electronic excitations (plasmons), and electron- 
hole (e-h) pairs. In metals, the excitation of e-h pairs 
leads to an excited or hot electron with energy above the 
Fermi level ep and to an excited or hot hole with en- 
ergy below ep- It is precisely the coupling of these hot 
carriers with the underlying substrate which governs the 
cross sections and branching ratios of electronically in- 
duced adsorbate reactions at metal surfaces. 

It is the intent of this Review to discuss the current 
status of the rapidly growing field of the dynamics of 
hot carriers in metals. The energy relaxation of these 
hot electrons and holes is almost exclusively attributed 
to the inelastic scattering with cold electrons below the 
Fermi level (e-e scattering) and with phonons (e-ph scat- 
tering), since radiative recombination of e-h pairs may 
be neglected. Assuming that the excess energy of the 
hot carrier is much larger than the thermal energy fc^T, 
the e-e scattering rate does not depend on temperature. 
Furthermore, for excitation energies larger than ^ 1 eV 
inelastic lifetimes are dominated by e-e scattering, e-ph 
interactions being in general of minor importance. Only 
at energies closer to the Fermi level, where the e-e inelas- 
tic lifetime increases rapidly, does e-ph scattering become 
important. 0,0 

Different techniques have recently become available for 
measuring hot-carrier lifetimes. Inverse photoemission 
(IPE) 6] and high resolution angle resolved photoemis- 
sion (ARPE)|2| provide an indirect access to the life- 
time of hot electrons and holes, respectively, by mea- 



suring the energetic broadening of transition lines after 
impigning an electron (IPE) or a photon (ARPE) into 
the solid. An alternative to IPE is two-photon photoe- 
mission (2PPE), 8J in which a first (pump) photon excites 
an electron from below the Fermi level to an intermediate 
state in the energy region e^? < e < £vac from where a sec- 
ond (probe) photon brings the electron to the final state 
above the vacuum level Evac- This technique can also be 
used to acces the lifetime of the intermediate state di- 
rectly in the time domain (time-resolved 2PPE),|^ by 
measuring the decrease of the signal as the probe pulse 
is delayed with respect to the pump pulse. Recently, 
it has been demonstrated that lifetime measurements of 
electrons and holes in metals can be done by exploiting 
the capabilities of scanniiig tunne ling microscopy (STM) 
and spectroscopy (STS),|l3, [H Il2| and ballistic elec- 
tron emission spectroscopy (BEES) has also shown to be 
capable of determi ning hot-electron relaxation times in 
solid materials. [13I Applications of these techniques 
include not only measurements of the scattering rates 
of hot carriers in solids, but also measurements of the 
lifetime of Shockley and image-potential states at metal 
surfaces. 

The early theoretical investigations of inelastic life- 
times and mean free paths of both low-energy electrons 
in bulk materials and image-potential states at metal 
surfaces have been described in previous reviews. jldllTf 
showing that they strongly depend on the details of the 
electronic band structure. [iJI Nevertheless, although ac- 
curate measurements of inelastic lifetimes date back to 
the mid 1990's, the first band-structure calculations of 
the e-e scattering in solids were not carried out until a 
few years ago.. 18, 19, 20, 21] 

Here we outline theoretical formulations of the hot- 
electron lifetime that is due to the inelastic scattering of 
the excited quasiparticle with occupied states below the 
Fermi level. Section 2 is devoted to the study of electron 
scattering processes in the framework of time-dependent 
perturbation theory. A discussion of the main factors 
that determine the decay of excited states is presented in 
Section 3, together with a review of the existing theoret- 
ical investigations of the lifetime of hot electrons in the 
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bulk of a variety of metals. The inclusion of exchange- 
correlation effects and chemical potential renormalization 
is addressed in Section 4, and the final Section presents 
an overview and future directions. 

Unless otherwise is stated, atomic units are used 
throughout, i.e., ~ h ^ rrie = 1. Hence, we use the 
Bohr radius, ao — /m^ = 0.529A, as the unit of length 
and the Hartree, H = e^/oo = 27.2 eV, as the unit of 
energy. The atomic unit of velocity is the Bohr veloc- 
ity, vq = ac = 2.19 X 10* cms~^, a and c being the fine 
structure constant and the velocity of light, respectively. 



THEORY 




We take a Fermi system of N interacting electrons at 
zero temperature (T = 0), and consider an external ex- 
cited electron interacting with the Fermi system. Fig. ^ 
depicts schematically a single inelastic scattering process 
for the excited hot electron. The hot electron in an ini- 
tial state 0i(r) of energy £i > e^? is scattered into the 
state 0/(r) of energy £/ {ep < £f < Si) hy exciting the 
cold Fermi system from its many-particle ground state of 
energy Eq to some many-particle excited state of energy 
En {En — Eq = Ei — Sf). By using the Fermi golden rule 
of time-dependent perturbation theory and keeping only 
the lowest-order term in the Coulomb interaction v{r, r') 
between the hot electron and the Fermi gas, the probabil- 
ity Pi^f per unit time for the occurrence of this process 
is found to bejj^ 

Pr->f = -2ldrl drV,*(r)0}(r') 

X Iml¥(r,r',£, -e/)0,(r')0/(r), (1) 
where W{r,r',uj) is the so-called screened interaction 

W{r,r';uj) = w(r,r') + J dridr2v(r, ri) 

X x(ri,i-2;w)w(r2,r') (2) 

and x(r, r'; uj) is the density- response function of the in- 
teracting Fermi system. ^JJ^j 

The total decay rate or reciprocal lifetime of the ex- 
ternal excited electron in the initial state 0i (r) of energy 
Ei is simply the sum of the probabilities Pi~,f over all 
available final states (j)f{r) with energies £/, i.e., 



(3) 



where the final states are subject to the condition ep < 

The single-particle wave functions and energies 0i j(r) 
and Eij entering Eq. |^ can be chosen to be the eigenfun- 
cions and eigenvalues of an effective Hartree, "231 Kohn- 
Sham,'23| or auasioarticle fisl l26| hamiltonian. 27] 




FIG. 1: Scattering of an external excited electron with a Fermi 
system of A*' interacting electrons at T = 0. The external 
electron in an initial state of energy Si > ef is scattered into 
an available state of energy e/ (ef < £/ < Ei) by exciting 
the cold Fermi system from its many-particle ground state of 
energy Eq to some many-particle excited state of energy E„ 

{En — Eo — Ei ~ Ef) 



Non-interacting Fermi sea 

If the Fermi sea is assumed to be a system of non- 
interacting electrons moving in an effective potential, in- 
stead of Eqs. Q-© one first considers the lowest-order 
probability P^^/ per unit time for the excited hot elec- 
tron in an initial state 4>i{r) of enery Ei to be scattered 
into the state 0/(r) of energy £/ by exciting one single 
electron of the Fermi sea from an initial state (j>i'{r) of 
energy £i' to a final state 0/' (r) of energy £/' (see Fig.©, 
and then obtains the probability P^/ by summing over 
all electrons in the Fermi sea (below the Fermi level) and 
all available final states for these electrons (above the 
Fermi level): 

^ 2 

Pi^f = 47r ^ n,/ (1 - Uf,) v{f, 6{Ei - £/ - £/- +Ei>), 

i'J' 

(4) 

where 



dr / dr'(b*{rm,{r')v{r,r')cj,f{r)cj,f,{r') (5) 



and Ui are Fermi-Dirac occupation factors, which at T 
are 



0{ef 



(6) 



9{x) being the Heaviside step function. The total decay 
rate is obtained by summing over all available final 
states of the excited hot electron, as in Eq. |j2Jl. 
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g of the Fermi sea can be described by the Thomas-Fermi 

^ ' approximation, 

O"* O iy^^(r,r') = f(r,r')e-«^^l""-'''l. (8) 

Here, qTF represents the Thomas- Fermi momentum, 
<1TF — (4'ZF/7r)^^^, qp being the Fermi momentum. 




FIG. 2: As in Fig. Q but now the Fermi system is assumed 
to be a system of non-interacting electrons. Hence, single- 
electron scattering events are represented, where the external 
excited electron of energy Si > sf is scattered into an avail- 
able state of energy e/ (ef < £f < £i), as in Fig. but now 
carrying one single electron of the Fermi sea from an initial 
state of energy e^/ (below the Fermi level) to a final state of 
energy Sfi (above the Fermi level). 



Alternatively, on can simply replace the imaginary part 
of the interacting density-response function x(r, r'; lo) en- 
tering Eq. (|2J) by its non-interacting counterpart: 

Im x°(r, r';Lu) = -2tt ^ n,, (1 - n},)<i>*, (r)0}, (r') 
■I' J' 

X(j)f,{v)(j),,{Y')5{uj-ef, +e^,). (7) 

Introduction of Eq. into Eqs. Q-Q yields again the 
probability of Eq. Q 

Due to the long-range of the bare Coulomb interaction 
?;(r, r') entering Eq. ||SJ), this approach yields a total de- 
cay rate that might be severely divergent, thereby re- 
sulting in a lifetime that would be equal to zero. How- 
ever, many-body interactions of the Fermi sea, which are 
fully included in the interacting density-response func- 
tion x(r, r'; 1^) entering Eq. are known to be responsi- 
ble for a dynamical screening of the Coulomb interaction 
leading to a finite lifetime of excited hot electrons. 

Many-body interactions of the Fermi sea are of- 
ten approximately introduced by simply replacing the 
long-range bare Coulomb interaction i;(r, r') entering 
Eq. |SJl by the frequency-dependent screened interaction 
of Eq. ((2Jl at the frequency uj = Ei — Sf, although this does 
not yield in general a result equivalent to that of Eq. 
In the long-wavelength limit, static (uj — > 0) screening 



Interacting Fermi sea: random-phase approximation 
(RPA) 

In an interacting Fermi sea, we need to compute the 
full interacting density-response function x(r, r';ct;) en- 
tering Eq. This function is also known to yield 
within linear-response theory the electron density in- 
duced in a many-electron system by an external potential 



(9) 



For many years, the dynamical screening in an inter- 
acting Fermi system has been succesfuUy described in a 
time-dependent Hartree or, equivalently, random-phase 
approximation (RPA). '2^ In this approach, the electron 
density induced by an external potential is obtained as 
the electron density induced in a non-interacting Fermi 
system by both the external potential y^^*(r, oj) and the 
potential l/*"''(r, oj) induced by the induced electron den- 
sity itself 



V""^{v,Lu)^ j dr'v{r,r')p'"'^{r',uj). 
Hence, in this approximation 

\r,uj)^ J dv' x°(r,r';L^) + j dn 



(10) 



And ( 



J rfi"2X°(i", ri; w)i)(ri, r2)x(r2, r'; uj) 



y-*(r',c.), 
(11) 

where r';w) is the density-response function of a 

system of non-interacting electrons. Comparing Eqs. © 
and one finds that in the RPA the interacting 

density-response function is obtained from the knowledge 
of the non-interacting density-response function by solv- 
ing the integral equation 



X(r,r';w) = r'; w) -f J dri y dra x"(i",i"i;w) 

X 'y(ri,r2)x(r2,r';w). (12) 



RESULTS 



General considerations 



The decay rate of excited electrons in an interacting 
Fermi sea is obtained from Eqs. the main ingre- 



dients being the hot-electron initial and final states [(f>i (r) 
and </'/(r)] and the imaginary part of the screened inter- 
action W{r,r';uj). lmW{r,r';uj) contains both a mea- 
sure of the probability for creating single-particle and 
collective excitations in a many-electron system and a 
measure of the screening of the interaction between the 
hot electron and the Fermi sea. In the case of low en- 
ergies, where collective modes cannot be produced, the 
hot-electron lifetime is mainly determined by a compe- 
tition between (i) the coupling of the initial state (f>i{r) 
with available states (f>f{r) above the Fermi level, (ii) the 
phase space available for the creation of electron-hole (e- 
h) pairs, and (iii) the dynamical screening of the Fermi 
sea. 
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FIG. 3: Schematic representation of the electronic band struc- 
ture of the (111) surface of Cu. 



Coupling with available states above the Fermi level 

The coupling of the excited electron (/)i(r) with avail- 
able states </>/(r) above the Fermi level strongly depends 
on whether the excited quasiparticle is a bulk or a surface 
state. 

A partially occupied band of Shockley surface states 
typically occurs in the gap of free-electron-like s,p 
bands. [iil Isof Another class of surface states occurs 
in the vacuum region of metal surfaces with a band gap 
near the vacuum level. These are unoccupied Rydberg- 
like image-potential states, which appear as a result of 
the self-interaction that an electron near the surface 
suffers from the polarization charge it induces at the 

surface, nil m 

Fig. 01 shows schematically the projection of the bulk 
band structure onto the (111) surface of the noble metal 
Cu. At the f point (A:|| = 0), the projected band gap 
extends from 0.9 eV below to 4.23 eV above the Fermi 
level, in such a way that both a Shockley (n = 0) and an 
image (n = 1) state are supported. 

Shockley surface states are known to be localized near 
the topmost atomic layer. However, image states are 
mainly localized outside the solid. Therefore, image 
states are expected to be weakly coupled with avail- 
able bulk states above the Fermi level and to live much 
longer than excited bulk states with the same energy: 
while the RPA broadening (or linewidth) of a bulk 
state at the energy of the n = 1 image state on Cu(lll) 
{si - Ef = 4.12 eV) is found to be 304 meV, the n = 1 
image-state linewidth is reduced to 29 eV. [231 

Coupling of the image state with the crystal occurs 
through the decay into bulk unoccupied states lying be- 
low the bottom of the projected band gap (which yields 
a linewidth of 17meV) and also through the decay into 
the unoccupied part of the n = Shockley surface state 
lying within the projected band gap (which leads to a 
linewidth of 12meV). A measure of the coupling of im- 
age states to bulk states of the solid can be given ap- 
proximately by the penetration of the image-state wave 



function into the solid, which in the case of the n = 1 
image state on Cu(lll) is found to be of 22%. However, 
the contribution to the linewidth coming from the decay 
into bulk states (17meV) is still well below 0.22 times 
the linewidth of a bulk state at the energy of the image 
state, i.e., 0.22 x 304 meV). This is due to the presence 
of a projected band gap at the surface, which in the case 
of Cu(lll) reduces considerably the phase space avail- 
able for real transitions of the n = 1 image state into 
bulk unoccupied states. This considerable reduction of 
the available phase space for real transitions does not oc- 
cur in the case of Cu(lOO); this explains the known fact 
that although the penetration of the n = 1 image state 
on Cu(lOO) is approximately 4 times smaller than in the 
case of Cu(lll) the ratio between the linewidths of the 
n ^ 1 image state on the (111) and (100) surfaces of Cu 
is smaller than 2. [3^ 



Phase space versus dynamical screening 

In the case of simple metals that can be described by 
a uniform free-electron gas of density no, which is char- 
acterized by the density parameter rg — {S/ATrnoY^^ /ao, 
a careful analysis of the phase space available for the 
creation of e-h pairs yields a decay rate of hot electrons 
with energies near the Fermi level (e^ — ep << £i) of the 
form [13 



Tr'^firs){e,-eFr. 

Furthermore, in the high-density limit (r^ 
® yieldHlIsllli 



(13) 

0), Eqs. 0- 



(14) 



This simple formula shows that the decay rate decreases 
as the electron density increases, which is also found to be 
true in the case of simple and noble metals (2 < < 6) 
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FIG. 4: Density of states of a free-electron gas with the 
electron density no equal to the average density of valence 
electrons in Al {rs = 2.07) and Na (r^ = 3.99). The 
RPA linewidth of high-density Al is found to be consider- 
ably smaller than the corresponding linewidth of low-density 
Na. 



and hot electrons with energies lying a few electronvolts 
above the Fermi level. 

As the electron density increases, the density of states 
(DOS) is known to increase. Hence, one might be 
tempted to conclude that as the electron density increases 
there are more electrons available for the creation of e-h 
pairs near the Fermi level, which would eventually yield 
an enhanced hot-electron decay. However, Eq. (|14|) shows 
that this is not the case. The reason for this behaviour 
is twofold. On the one hand, due to momentum and 
energy conservation the number of states available for 
real transitions in metals is typically weakly dependent 
on the actual electron density. On the other hand, as 
the electron density increases the ability of electrons to 
screen the Coulomb interaction with the external excited 
electron also increases, which leads to a smaller screened 
interaction and a reduced hot-electron decay. 

Fig. 21 shows the DOS of a free-electron gas with the 
electron density tiq equal to the average density of va- 
lence electrons in the simple metals Al (rj = 2.07) and 
Na {vs — 3.99), together with the corresponding RPA 
linewidths [as obtained from Eqs. Q-©] of hot electrons 
with energy lying 1 eV above the Fermi level. While in 
the case of Al, with a large DOS near the Fermi level, the 
linewidth is 14 meV, the linewidth of hot electrons with 
£i — ep = 1 eV in Na is 59meV. Hence, hot electrons 
live longer in Al than in Na, which is basically due to 
the strong screening characteristic of high-density metals 
like Al. 



Free-electron gas 

For many years, theoretical predictions of the electron 
dynamics of bulk states in solids had been based on a free- 
electron gas (FEG) or jellium description of the solid, in 
which a homogeneous assembly of interacting electrons 
is assumed to be immersed in a uniform positive back- 
ground. In this model there is translational invariance, 
the one-particle states entering Eq. Q are momentum 
eigenfunctions, and Eqs. CJ-© are easily found to yield 

^k' = -2/(03l^^<i.- (15) 

where the energy transfer to = ek~ek-q [here, = 
is subject to the condition 0<a;<ek — £_f, and 

Vq and Xq,u/ being Fourier transforms of the bare 
Coulomb interaction v{r,r') and the interacting density- 
response function x(r, r';^), respectively. The screened 
interaction Wq^ai is usually expressed in terms of the in- 
verse dielectric function e"^^, as follows 

Wq.^^Vqe-1,, (17) 

where 

eq,L = 1 + Xq,t^ Vq. (18) 

Non-mteracting Fermt sea 

If the Fermi sea is assumed to be a system of non- 
interacting electrons, Xslu reduces to the well-known 
Lindhard function Xq,wjlH)ll3- One can then write the 
imaginary part of the screened interaction of Eq. (|16|l as 

/dk' 
j—^ni,>{l-nk'+q) 

X |Wq|^(5(cj -tJk'+q+^k')- (19) 

Interacting Fermi sea: RPA 

In the RPA, the interacting density-response function 
Xq.uj is obtained from Eq. ifT^ . i.e., 

(20) 

or, equivalently [see Eq. ((TH|l ]. 

eq,(^ = 1 -Xq,c:.^'q- (21) 

Introducing Eq. H20() into Eq. H16|l . one finds 

/dk' 
^^^^3 »-k'(l - "k'+q) 

X |Wq,^|^5(w -tJk'+q +Wk'), (22) 
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FIG. 5: RPA lifetime of hot electrons in a free-electron 
gas, divided by tqf [Eqs. 1)13^ and l|14|l ]. as a function of 
the electron-density parameter and for hot electrons in the 
vicinity of the Fermi surface. 



FIG. 6: RPA lifetime of hot electrons in a free-electron gas, 
divided by tqf [Eqs. 1131 and 114ll l. as a function of et — £f 
for an electron density equal to that of valence electrons in Al 
{rs = 2.07) 



which is of the form of Eq. (|19|l with the bare Coulomb 
interaction Vq replaced by the RPA screened interaction 



W, 



Beyond the RPA, ImWq.i^ cannot always be ex- 



pressed as in Eq. 

We note that in the high-density limit (r^ — > 0) and 
for hot electrons with energies lying near the Fermi level 
(Ck - Ef << ^k), introduction of Eq. into Eq. ifT^ 
yields the decay rate given by Eqs. and (|14|l . which 
will be referred as Tgp. 

Figs. 6 and 7 show the RPA lifetime Tk of hot electrons 
in a free-electron gas, divided by tqp, as a function of 
the electron-density parameter for hot electrons in the 
vicinity of the Fermi surface (Fig. [SJ , and as a function 
of Ek — £_F for an electron density equal to that of va- 
lence electrons in Al (Fig. 0). Although the high-density 
limit of Eq. H14() only reproduces the full RPA calcula- 
tion as rs — > 0, Fig. shows that differences between 
this high-density approximation (tqp) and the full RPA 
calculation (r) are very small at electron densities with 

< 2 and go up to no more than 7% at = 6. In the 
inset of Fig. [5l the RPA lifetime is represented, showing 
that as occurs in the high-density limit [see Eq. H14|l ] the 
RPA lifetime Tk increases very rapidly with the electron 
density. 

In the limit £k — > £_f the available phase space for real 
transitions is simply Ek — £_f, which yields the (sk ~ ef)^ 
quadratic scaling of Eq. l(T^ . However, as the energy 
increases momentum and energy conservation prevents 
the available phase space from being as large as Ek — 
ep- As a result, the actual lifetime departures from the 
quadratic scaling predicted for electrons in the vicinity 
of the Fermi level (see Fig. EJ, differences between the 
full RPA lifetime Tk and the lifetime tqp dictated by 
Eqs. H13|l and 114|l ranging from « 0.5% at Sk ~ £_f to 
w 40% at ek -£f = 4eV. 



Random-fc approximation 

For the description of lifetimes and inverse mean free 
paths (IMFP) in non-free-electron solids, several au- 
thors have employed the so-called random-fc approxi- 
mation first considered by Berglund and Spicer'38] and 
by Kane.'39l[ The starting point of this approximation 
is Eq. Q with the bare Coulomb interaction v{r,r') 
replaced by a frequency-dependent screened interac- 
tion iy(r,r';w). The random-A; approximation is then 
the result of replacing all the squared matrix elements 
i'->f' ^ 

^iJ^f (uj) by their average Mi over all available states 

/', and /. All the summations entering Eqs. (O 
and can then be replaced by integrations over the 
corresponding density of states, and one can write the 
reciprocal lifetime as 



eF-(ei-e) 



M, / dep{e) I ds' p{£')p{e^+e' -e), 

(23) 

where the quantities p(£) are assumed to include spin, 
i.e., 



d£p{e) = N, 



(24) 



N being the total number of electrons. 

A further simplification may be achieved if the density 
of states is assumed to take the constant value pi below 
and above the Fermi level. Eq. (|23() then reduces to 



■ pf {e, - ep)' 



(25) 



As the electron density increases, both the DOS and the 
screening of the Fermi sea increase. Since an enhanced 
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screening yields a reduced Mi factor, the actual depen- 
dence of the hot-electron lifetime on the electron den- 
sity would be the result of the competition between DOS 
and screening effects. Screening effects typically domi- 
nate leading to a lifetime that increases with the electron 
density, as occurs in the case of a free-electron gas. 

The random-fc approximation was first used by 
Berglund and Spicer in order to explain experimental 
photoemission studies of Cu and Ag.[2S| Scattering rates 
of electrons in Si were computed by Kane.js^ and Kro- 
likowsky and Spicer |4oj employed the random-fc approx- 
imation to calculate the energy dependence of the IMFP 
of electrons in Cu from the knowledge of density-of-state- 
distributions in this material which had been deduced 
from photoelectron energy-distribution measurements. 

More recently, Penn et al. used the random-fc approxi- 
mation to analyze the spin-polarized electron-energy-loss 
spectra and hot-electron lifetimes in ferromagnetic Fe, 
Ni, Co, and Fc-B-Si alloys. [ill li^ The experimental spin- 
dependent hfetime of the n = 1 image state on Fe(llO) 
was also interpreted by using this approximation. |43 
Drouhin developed a model to evaluate the scattering 
cross section and spin-dependent inelastic mean free path 
from the knowledge of density-of-state distributions, and 
applied it to Cr, Fe. Co. Ni, Gd, Ta, and the noble metals 
Cu, Ag, and Au.|4J| An analogous approach was devel- 
oped by Zarate et al., where simple approximations to 
the DOS were used to obtain analytical expressions for 
the electron lifetimes in transition metals. [45j 

The random-fc approximation has also been discussed 
and compared to first-principles calculations by Zhukov 
et a/..[4y, |47| It was shown that when initial and final 
states are either both sp or both d states then with 
a proper choice of the averaged matrix elements Mi 
Eq. (|23|l yields lifetimes in reasonable agreement with 
more elaborated calculations. 



First-principles calculations 

First-principles calculations of the scattering rates of 
hot-electrons in periodic solids were first carried out only 
a few years ago by Campillo et In this work, hot- 

electron inverse lifetimes were obtained from the knowl- 
edge of the on-shell electron self-energy of many-body 
theory, which in the so-called G'^W approximation yields 
exactly the same result as Eqs. O"® above. 

For periodic crystals, the single-particle wave functions 
entering Eq. are Bloch states </'k,i(r) and 0k-q./(r) 
with energies £k,i and Ek-q,/, i and / representing band 
indices. Hence, introducing Eq. ^ into Eq. (0) and 
Fourier transforming one finds 

f ^ ' G.G' 

X B,j(q+G')ImI^cj,G'(q,w), (26) 



where the energy transfer oj = ek,i — £k-qj is subject 
to the condition < a; < £k,i — ef, the integration is 
extended over the first Brillouin zone (BZ), the vectors 
G and G' are reciprocal lattice vectors, Bif represent 
matrix elements of the form 



i3,/(q + G) = j dr0*,,(r)e'('i+«)--0k-q,/( 



(27) 



and Wg^g' (q, u) are Fourier coefficients of the screened 
interaction W{y, r'; w). As in the case of the free-electron 
gas, the Fourier coefficients Wg.g' (q, w) are usually ex- 
pressed in terms of the inverse dielectric matrix 



W^G,G'(q,t^) ='y(q + G)eQ|g,(q,tj), 



(28) 



ti(q) being the Fourier transform of the bare Coulomb 
interaction u(r,r'). In the RPA, 

ec.G' (q, oj) = 5g^G' - Xg,G' (q, ^) «(q + G'), (29) 

where x%, g' (q^ ^) ^'^'^ Fourier coefficients of the non- 
interacting density-response function x'^(r,r';w). 

Couplings of the wave vector q + G to wave vectors 
q -I- G' with G ^ G' appear as a consequence of the 
existence of electron-density variations in the solid. If 
these terms, representing the so-called crystalline local- 
field effects, are neglected, one can write Eq. (|26|) as 

dq 



'k,i 



BIM + G) 



X B,j(q + G)i;(q + G 



Im£G,G(q,^) 
|eG,G(q,t^)|^ 



(30) 



This expression accounts explicitly for the three main in- 
gredients entering the hot-electron decay process. First 
of all, the coupling of the hot electron with available 
states above the Fermi level is dictated by the matrix 
elements Bij{q + G). Secondly, the imaginary part of 
the dielectric matrix eG,G(q, w) represents a measure of 
the number of states available for the creation of e-h 
pairs with momentum and energy q -I- G and u>, respec- 
tively. Thirdly, the dielectric matrix in the denominator 
accounts for the many-body e-e interactions in the Fermi 
sea, which dynamically screen the interaction with the 
external hot electron. 

Since for a given hot-electron energy e there are in 
general various possible wave vectors and bands, one can 
define an energy-dependent reciprocal lifetime by doing 
an average over wave vectors and bands. As a result of 
the symmetry of Bloch states, one finds t^^- = t^], with 
S representing a point group symmetry operation in the 
periodic crystal. Thus, one can write 



IBZ 



'k,n 



n Lk "^k.n 



(31) 



where mk,n represents the number of wave vectors k lying 
in the irreducible element of the Brillouin zone (IBZ) with 
the same energy e. 
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FIG. 7: Hot-electron lifetimes in Al, as a function of the hot- 
electron energy e with respect to the Fermi level ef- Solid 
circles represent the ab initio RPA calculation, as obtained 
after averaging the lifetime broadening of Eq. (12611 over 
all states with the same energy e. The solid line represents 
the RPA lifetime of hot electrons in a FEG with Ts = 2.07. 
Open circles represent the result obtained from Eq. (1261 by 
replacing the hot-electron initial and final states entering the 
coefficients Bif by plane waves, but with full inclusion of the 
band structure in the evaluation of ImVKG^c (q, uj). The inset 
exhibits scaled lifetimes of hot electrons in Al. Solid circles 
and the solid line represent first-principles and FEG calcula- 
tions, respectively, both within RPA. The dashed line repre- 
sents the prediction of Eqs. 11311 and 1141 . 

Plane-wave (PW) basis 

The calculations reported in Refs. 0, l^^, 0| for the 

lifetime of hot electrons in the simple metals Al, Mg, and 
Be, and the noble metals Cu and Au were carried out 
from Eqs. H2t)|) - (|29|l by expanding all one-electron Bloch 
states in a plane-wave (PW) basis: 

</'M(r) = ^ ^^ik,,:(G)e'(''+°)•^ (32) 

G 

where fl represents the normalization volume. More re- 
cently, similar calculations were reported by Bacelar et 
al. for the lifetime of hot electrons in six transition met- 
als: two fee metals (Rh and Pd), two bee metals (Nb and 
Mo), and two hep metals (Y and Ru).|4^ 

In this approach, one first solves self-consistently 
for the coefficients Mk,i(G) the Kohn-Sham equa- 
tion of density-functional theory (DFT), with use of 
the local-density approximation (LDA) for exchange 
and correlation l49i a nd non-local norm-conserving ionic 
pseudopotentialsjS^ to describe the electron-ion inter- 
action. From the knowledge of the eigenfunctions and 
eigenvalues of the Kohn-Sham hamiltonian, one can pro- 
ceed to evaluate the non-interactin g de nsity-response ma- 
trix Xg G'i^y^) (see, e.g., Ref. [I3) and the dielec- 
tric matrix eG.G'(q, w), a matrix equation must then be 



solved for the inverse dielectric matrix, and the scattering 
rate is finally computed from Eq. (|26() with full inclusion 
of crystalline local-field effects. 

Here, we are showing the results of ab initio RPA 
calculations of the average lifetime of hot electrons in 
the simple metal Al and the noble metals Cu and Au, 
all obtained from Eqs. and (|^ . and we com- 

pare these results to the RPA lifetime of hot electrons 
in the corresponding FEG as obtained from Eqs. (jl5|l 
and (jl6(l . Both ab initio and FEG calculations were per- 
formed within the very same many-body framework, and 
the comparison between them indicates that while band- 
structure effects reduce the lifetime of hot electrons in 
Al, the presence of non-free-electron-like d bands in the 
noble metals considerably enhances the lifetime of hot 
electrons in Cu and Au. 

In Fig. [7| we show ab initio (solid circles) and FEG 
(solid line) RPA calculations of the average lifetime r(e) 
[see Eq. (|31|) ] of hot electrons in the face-centered-cubic 
(fee) Al, as a function of the hot-electron energy e with 
respect to the Fermi level e f ■ These calculations show 
that even in a free-electron metal like Al band-structure 
effects play a key role lowering the hot-electron lifetime 
by a factor of ~ 0.65 for all electron energies under study. 
In order to understand the origin of band-structure ef- 
fects, an additional calculation is represented by open cir- 
cles where the hot-electron initial and final Bloch states 
entering the coefficients Bij of Eq. H27I) have been re- 
placed by plane waves (as in a FEG) but keeping the full 
inverse dielectric matrix of the crystal. This calculation, 
which lies nearly on top of the FEG curve, shows that 
band-structure effects on both e-h pair creation and the 
dynamical screening of the Fermi sea are very small, as 
occurs in the case of slow ions.lsj However, the coupling 
of hot-electron initial and final Bloch states appears to 
be very sensitive to the band structure, which in the case 
of Al shows a characteristic splitting over the Fermi level 
thereby opening new channels for electron decay and re- 
ducing the lifetime. [20j 

Similar results to those exhibited in Fig. |7|for Al were 
obtained for the hexagonal closed-packed (hep) Mg.[2(il| 
whose band structure also splits just above the Fermi 
level along certain symmetry directions. However, the 
splitting of the band structure of this material is not as 
pronounced as in the case of Al, and the departure of the 
hot-electron lifetime in Mg from the corresponding FEG 
calculation with Ts = 2.66 was found to be of about 25%, 
smaller than in Al. 

Figs. 181101 exhibit ab initio (solid circles) and FEG 
(solid lines) RPA calculations of the average lifetime r(e) 
of hot electrons in the fee metals Cu and Au, again as a 
function of the hot-electron energy e with respect to the 
Fermi level. Cu and Au are noble metals with entirely 
filled 3d and 5d bands, respectively. Slightly below the 
Fermi level, a,t e — ep ^ '2 eV, we have d bands capable 
of holding 10 electrons per atom, the one remaining elec- 




FIG. 8: Hot-electron lifetimes in Cu, as a function of the hot- 
electron energy e with respect to the Fermi level ef- Solid 
circles represent the ab initio RPA calculation, as obtained 
after averaging the lifetime broadening of Eq. (12611 over 
all states with the same energy e. The solid line represents the 
RPA lifetime of hot electrons in a FEG with =2.67. The 
dotted line represents the ab initio RPA calculation of the av- 
erage lifetime r(£), but with the 3d shell assigned to the core 
in the pseudopotential generation. The squares represent the 
result of replacing in Eq. 13011 all one-electron Bloch states by 
plane waves (FEG calculation) but keeping the actual density 
of states in the evaluation of Imec c' {^., Open circles rep- 
resent the result obtained from Eq. 1301 by replacing the hot- 
electron initial and final states enetring the coefficients Bif 
by plane waves and the dielectric function in |eG,G(q, 
by that of a FEG with Ts — 2.67, but with full inclusion of 
the band structure in the calculation of Imec^G' (q, '^)- Tri- 
angles represent the result obtained from Eq. 13011 by replac- 
ing the hot-electron initial and final states entering the co- 
efficients Bif by plane waves, but with full inclusion of the 
band structure in the evaluation of both Imec (q, tj) and 
|eG,G(q,c^)r^. 



tron being in both Cu and Au in a free-electron-like band 
below and above the d bands. Hence, a combined descrip- 
tion of both delocalized s valence bands and localized d 
bands is needed to address the actual electronic resi; 
of these metals. The results reported in Refs. 
and presented in Figs. I8I10I where found by keeping all 
4s ^ and 3(i^° Bloch states (in the case of Cu) and all 
6s ^ and Sd^" Bloch states (in the case of Au) as valence 
electrons in the generation of the pseudopotential. 

Also represented in Fig. |S| are additional calculations 
for the hot-electron lifetime in Cu, which help to under- 
stand the origin and impact of band-structure effects in 
this material. These are an ab initio RPA calculation 
similar to the full RPA calculation represented by solid 
circles, but with the 3d shell assigned to the core (dashed 
line), and three approximated calculations from Eq. H3l)|l 
[thus neglecting crystalline local-field effects] in which the 
actual band structure is only considered in (i) the DOS 
entering the evaluation of Imecclq, w) (squares), (ii) 
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FIG. 9: Experimental lifetimes of hot electrons in Cu, multi- 
plied by (e — gf)^ and as a function of (e — £f). The open cir- 
cles represent the TR-2PPE measurements reported by Knoe- 
sel et al. \QQ\ for the lifetime of very- low-energy hot electrons 
in Cu(lll). The TR-2PPE measurements reported by Ogawa 
et al. |52l for the lifetime of hot electrons in Cu(llO) are repre- 
sented by open squares. For comparison, the ab initio calcula- 
tions represented in Fig.|H]by solid circles are also represented 
in this figure (solid circles), together with the lifetimes of hot 
electrons in a FEG with Ts — 2.67 (solid line) and the corre- 
sponding approximation of Eqs. 1131 and 1141 dashed line). 



the full evaluation of ImeG,G(q, ^) (open circles), which 
also contains the coupling between states below and 
above the Fermi level entering the production of c-h pairs, 
and (iii) the full evaluation of the imaginary part of the 
inverse dielectric function ImeG,G(qi'^)/ |eG,G(q, w)|^ 
(triangles). 

An inspection of Fig.|51shows that when the 3d shell is 
assigned to the core ab initio calculations (dashed line) 
nearly coincide with the FEG prediction (solid line). This 
must be a consequence of the fact that band-structure ef- 
fects in Cu are nearly entirely due to the presence of lo- 
calized d electrons. In the presence of d electrons, there 
are obviously more states available for the creation of 
e-h pairs, and one might be tempted to conclude that 
d electrons should yield an enhanced hot-electron decay, 
especially at the opening of the c?-band scattering channel 
at about 2 eV below the Fermi level. Indeed, this is pre- 
cisely the result of calculation (i) represented by squares, 
which is very close to the calculation reported by Ogawa 
et aL. [53,153 Nevertheless, a full ab initio evaluation of 
ImeG,G(q, 1^) yields the calculation (ii) represented by 
open circles; this shows that there is no coupling between 
d and sp electrons below and above the Fermi level and 
there is, therefore, little impact of d electrons on the pro- 
duction of e-h pairs. [H^l The key role that d electrons play 
in the hot-electron decay is mainly due to screening ef- 
fects. The presence of d electrons gives rise to additional 
screening, thus increasing the lifetime of hot electrons for 
all excitation energies, as shown by the result of calcu- 
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FIG. 10: Hot-electron lifetimes in An, as a function of the hot- 
electron energy e with respect to the Fermi level ef- Solid 
circles represent the ab initio RPA calculation, as obtained 



in Ref. 



after averaging the lifetime broadening ^ of 



Eq. 12611 over all states with the same energy e. The open 
squares represent the experimental measurements of Ref. |6^. 
The solid line represents the RPA lifetime of hot electrons in 
a FEG with rs — 3.01. The triangles represent the result 
obtained from Eq. 1301 by replacing the hot-electron initial 
and final states entering the coefficients Bif by plane waves 
and the dielectric function in |eG,G(q,a;)|~^ by that of a FEG 
with Vs — 3.01, but with full inclusion of the band structure 
in the calculation of Imec (q, tj). The inset exhibits scaled 
lifetimes of hot electrons in Au. The solid circles and the 
solid line represent first-principles and FEG calculations, re- 
spectively, both within RPA. The dashed line represents the 
prediction of Eqs. 11311 and 1141 . The triangles represent the 
result obtained from Eq. 13011 by replacing the hot-electron 
initial and final states entering the coefficients Bif by plane 
waves and the dielectric function in |eG,G(q,a;)|~^ by that of 
a FEG with Vs = 3.01, but with full inclusion of the band 
structure in the calculation of ImeG^Q/ (q, oj). 



lation (iii) represented by triangles which includes the 
screening of d electrons in the ab initio evaluation of the 
denominator of Eq. (|30|l . The screening of d electrons is 
found to increase hot-lectron lifetimes by a factor of 3. 
Finally, differences between calculation (iii) (triangles) 
and the full ab initio calculation (solid circles), which 
are only visible at hot-electron excitation energies very 
near the Fermi level, are due to a combination of band- 
structure effects on the hot-electron initial and final sp 
states above the Fermi level and crystalline local-field ef- 
fects not included in Eq. I|30|l . 

The first time-resolved 2PPE (TR-2PPE) experiments 
on Cu were performed by Schmuttenmaer et aL.js^ The 
electron dynamics on copper surfaces was later investi- 
gated by several groups.js^jEEIUIililiSE^ Although 
there were some discrepancies among the results obtained 
in different laboratories, most measured lifetimes were 
found to be much longer than those of hot electrons in 
a free gas of 4s^ valence electrons. The first-principles 



(111), 
in 



calculations reported in Ref. 0| (see also Fig. ^ indi- 
cate that this is mainly the result of the screening of d 
electrons. 

Fig. |51 exhibits (open circles) the TR-2PPE measure- 
ments reported by Knoesel et aL[60|| for the lifetime of 
very- low-energy electrons in Cu(lll). Since the energy of 
these electrons is less than ~ 2eV above the Fermi level, 
d electrons (the d-band threshold is located at ~ 2 eV be- 
low the Fermi level) cannot participate in the creation of 
electron-hole pairs but can participate in the screening of 
e-e interactions, thereby increasing the hot-electron life- 
time. Fig. |51 shows that the measurements of Knoesel et 
al. are in reasonable agreement with the first-principles 
calculations reported in Ref. 0| (solid circles). The re- 
laxation dynamics at the low-index copper surfaces 
(100), and (100) were investigated by Ogawa et al. 
a wider energy range. Since the experiment detects the 
photoemitted electrons in the normal direction (ky — 0), 
we have plotted in Fig. (open squares) the measured 
lifetimes of hot electrons at the (110) surface of Cu, the 
only surface with no band gap for electrons emitted in 
this direction. At large electron energies (e — > 2 eV) , 
there is very good agreement between calculated lifetimes 
(solid circles) and the measurements reported in Ref. [s^ 
(open squares). Nevertheless, the calculations cannot ac- 
count for the large increase in the measured lifetimes re- 
ported in Ref. |52l | at lower electron energies. The origin 
of this discrepancy was discussed by Schone et aZ.|6ll|. 
who concluded that the increase in the experimentally 
determined lifetime might be related to the presence of 
a photohole below the Fermi level leading to transient 
excitonic states. 

Ab initio RPA average lifetimes in Au, as reported 
in Ref. j23|, are represented in Fig. ^| (solid circles), 
together with accurate TR-2PPE measurements (open 
squares) where the relaxation from electron transport to 
the surface was expected to be negligible, [g^ This figure 
shows that there is very good agreement between theory 
and experiment, which are both found to be nearly five 
times larger than the lifetime of hot electrons in a free gas 
of 6s^ valence electrons (solid line). Nevertheless, more 
accurate linearized augmented plane- wave (LAPW)|Q| 
and plane- wave calculations of the lifetime of hot elec- 
trons in Au have been carried out recently, which yield 
smaller lifetimes than those reported in Ref. [ij and rep- 
resented in Fig. ^1 by an overall factor of ~ 1.4. These 
calculations have been found to accurately reproduce the 
BEES spectra for the two prototypical Au/Si and Pd/Si 
sy stems lUj. 

The triangles of Fig. ^| show the result obtained in 
Ref. [iJl from Eq. \60\ (thus neglecting crystalline local- 
field effects) by including the actual band structure of Au 
in the evaluation of Im eG_G(q, cj) but treating the coef- 
ficients Bij- and the denominator |eG.G(q, 'j-')|^ as in the 
case of a FEG with rg = 3.01. The result of this calcu- 
lation is very close to the FEG calculation (solid lines), 



which shows that the opening of the d-band scattering 
channel for e-h pair production does not play a role. As 
in the case of Cu, differences between the FEG and the 
full ab initio calculation are mainly a consequence of vir- 
tual interband transitions due to the presence of d elec- 
trons, which give rise to additional screening and largely 
increase the hot-electron lifetime in Au. 

The change 6e in the real part of the long-wavelength 
dielectric function that is due to the presence of d elec- 
trons in the noble metals is known to be practically con- 
stant at the low frequencies involved in the decay of low- 
energy hot electrons in these materials . |63l Im . l65l| Hence, 
in order to investigate hot-electron lifetimes and mean 
free paths in the noble metals Quinnlo^ treated the FEG 
as if it were embedded in a medium of dielectric constant 
eo = l+Se instead of unity. The corrected lifetime is then 

1/2 

found to be larger by a factor of , i.e., ~ 2.5 for both 
Cu and Au. A more accurate analysis has been carried 
out recently, in which the actual frequency dependence 
of eo is included,!^ showing that for the low energies 
of interest the corrected lifetime is still expected to be 
enhanced by a factor of ~ 2.5 for both Cu and Au. Nev- 
ertheless, the ab initio calculations exhibited in Figs. 9-11 
indicate that the role that occupied d states play in the 
screening of e-e interactions is much more important in 
Au than in Cu. 

Finally, we note that the screening of d electrons does 
not depend on whether the hot electron can excite d elec- 
trons (both in Cu and Au the d-band scattering channel 
opens at ~ 2eV below the Fermi level) or not. Thus, 
d-screening effects do not depend on the hot-electron en- 
ergy, and average lifetimes are therefore expected to ap- 
proximately scale as (e — as in the case of a FEG. 
Nonetheless, departures from this scaling behaviour can 
occur when the hot-electron lifetimes are measured along 
certain symmetry directions, as discussed in Ref. |l^ 



Linear mujfin-tin orbital (LMTO) basis 

For noble and transition metals containing d and / 
electrons it is sometimes more convenient to perform 
full all-electron calculations based on the use of local- 
ized basis like linear mufhn-tin orbital (LMT0),p3| lin- 
ear combination of atomic orbital (LCAO),!^^ or LAPW 
basis. [l^ The lifetime of hot electrons in the noble met- 
als Cu, Ag, and Au, and the Ad transition metals Nb, 
Mo, Rh, and Pd was determined by Zhukov et aL0i0| 
by using an LMTO basis in the so-called atomic-sphere 
approximation (ASA), in which the Wigner-Seitz (WS) 
cells are replaced by overlapping spheres. 

In the LMTO method, one-electron Bloch states are 
expanded as follows 
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FIG. 11: Hot-electron lifetimes in Nb and Pd, as a function of 
the hot-electron energy e with respect to the Fermi level ef- 
Solid and open circles represent ab initio RPA calculations 
for Nb and Pd, respectively, as obtained after averaging the 
lifetime broadening of Eq. 12611 over all states with the 
same energy e. 



where L — (Z,m) represents the angular momentum and 
XR,L(r, k) are the LMTO basis wave functions, which 
in the ASA are constructed from the solutions </3R„;(r) 
to the radial Schrodinger equation inside the overlap- 
ping muffin-tin sphere at site R and their energy deriva- 
tives. Hot-electron reciprocal lifetimes can then be evalu- 
ated from Eqs. (|26|I - H29|I by first solving self-consistently 
for the coefficients Uk. j(R', ^) the Kohn-Sham equation 
of DFT with the use of the LDA for exchange and 
correlation 49j and then computing the dielectric matrix 
^G,G'{^.,^) from first principles. 

In Fig. ^] we plot the ah initio RPA calculations that 
we have carried out from Eqs. for the average 

lifetime T(e) [see Eq. H31|l ] of hot electrons in the body- 
centered-cubic (bcc) Nb and the fee Pd. Nb (4d'*5s^) 
is a transition metal with a partially filled 'id band and 
one 5s valence electron per atom. In this material, the 
bands of d states lie in the energy interval from 4eV 
below to 6eV above the Fermi level. Hence, contrary 
to the case of the noble metals where hot electrons de- 
cay by promoting sp electrons from below to above the 
Fermi level, in the case of Nb and other transition metals 
d electrons also participate in the creation of e-h pairs. 
Accordingly, the net impact of d electrons in the decay 
of hot electrons in Nb will be the result of the competi- 
tion between the opening of a new scattering channel and 
the presence of d-electron screening. Valence (5s^) elec- 
trons in Nb form a FEG with electron-density parameter 
Ts = 3.07. However, the actual lifetimes of hot electrons 
in Nb are considerably shorter than those of hot electrons 
in the free gas of 5s^ valence electrons, which is obviously 
due to the presence of d electrons. Furthermore, if one 
naively considered a FEG with 5 valence electrons per 
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atom {vs = 1.80), one would predict a lifetime that is 
- (3.07/1.80)^/2 4 times larger [see Eq. (O] than in 
the case of a FEG with = 3.07 and nearly 8 times 
larger than predicted by our ab initio calculation. On 
the one hand, d electrons below and above the Fermi level 
in real Nb strongly couple for the creation of e-h pairs, 
thereby strongly decreasing the hot-electron lifetime for 
all electron energies under consideration. On the other 
hand, screening of localized d electrons is much weaker 
than in the case of free electrons. Hence, while screening 
effects dominate in the case of a FEG, the effect of the 
DOS near the Fermi level dominates in the case of a gas 
of electrons with a large number of d states, which yields 
the short average lifetimes plotted in Fig. ^2 

Pd (4(i^") is a transition metal with a filled Ad band 
and no valence sp electrons. In this material, the bands of 
d states mainly lie below the Fermi level, though states 
above the Fermi level still have a small but significant 
d component. The average density of 10 electrons per 
atom in Pd is higher than that of 5 electrons per atom 
in Nb; a FEG picture of the hot-electron dynamics in 
these materials would, therefore, lead to longer lifetimes 
in Pd. Nevertheless, our ab initio calculations show that 
in the case of low-energy electrons with energies below 
4eV this is not the case, i.e., hot electrons in Pd live 
shoter than in Nb. This is again due to the fact that as 
far as d electrons are concerned the effect of the density 
of d states near the Fermi level shortening the lifetime is 
more important than the effect of c?-screening increasing 
the lifetime, simply because of the limited ability of d 
electrons to screen the e-e interactions. The density of d 
states near the Fermi level is higher in Pd than in Nb, 
thereby leading to lifetimes of low-energy electrons that 
are shorter in Pd than in Nb. At energies e — ep > 4eV 
e-h pair production is restricted to the lowest energies, 
due to the fact that for the highest energies d electrons 
below the Fermi level do not couple with sp electrons at 
a few eVs above the Fermi level. Thus, at these energies 
e-h pair production in Pd is not as efficient as in the case 
of Nb, which yields hot-electron lifetimes that are larger 
in Pd than in Nb. 



GREEN FUNCTION FORMALISM 

In the theoretical framework of section II we have re- 
stricted our analysis to the case of an external excited 
electron interacting with a Fermi system of interacting 
electrons. Nevertheless, a more realistic analysis should 
include the excited hot electron as part of one single 
Fermi system of A" -I- 1 interacting electrons in which an 
electron has been added in the one-particle state 4>i (r) at 
time t' . In the framework of Green function theory, [23 
the probability that an electron will be found in the same 
one-particle state at time t > t' is obtained from the 
knowledge of the one-particle Green function of a system 



of A' interacting electrons. One finds. [2 6| 

P,it',t) ^exp[2ImE,it~t')], (34) 

Si being the so-called quasiparticle energy, i.e., the pole 
of the Green function. Hence, the total decay rate or 
reciprocal lifetime of the quasiparticle is simply 



~2Im£', 



(35) 



The quasiparticle energy Ei can be approximately ex- 
pressed in terms of the electron self-energy [I](r, r'; £";)] 
and the eigenfunctions [(/'i(r)] and eigenvalues [si] of 
an effective single-particle hamiltonian [H = — Vr/2 + 
V{r)]-M 



(36) 



where 



AE(£;,) = J dv J dr'^*{r)[j:{r,r';E,) 

- V{r)6ir^r')]Mr'). (37) 

The total decay rate or reciprocal lifetime of the quasi- 
particle can, therefore, be obtained as follows 

rri = -2ImAS(S,) 

= -2 J drj drV*(r)ImE(r,r';i;O0.(r').(38) 

Within many-body perturbation theory,^] it is pos- 
sible to obtain the electron self-energy T,{r,r';Ei) as a 
series in the bare Coulomb interaction w(r,r'), but due 
to the long range of this interaction such a perturba- 
tion series contains divergent contributions. Therefore, 
the electron self-energy is usually rewritten as a series 
in the time-ordered screened interaction W'^^{r,r';u!), 
which can be built from the knowledge of the Green func- 
tion itself and coincides for positive frequencies (w > 0) 
with the retarded screened interaction of Eq. In 
the RPA, the Green function entering the screened in- 
teraction W'^'-^{r,r';uj) is replaced by its noninteracting 
counterpart and for positive frequencies one simply finds 
Eq. with the density-response function of Eq. H12|l . 

To lowest order in the screened interaction, the self- 
energy E(r, r'; Ei) is obtained from the product of the 
Green function and the screened interaction, and is there- 
fore called the GW self-energy. If one further replaces 
both the Green function G(r, r';uj) and the quasiparticle 
energy Ei entering the GW self-energy by their nonin- 
teracting counterparts [AT,{E,) AS°(ei)], Eq. ^ 
is easily found to yield an expression for the decay rate 
which exactly coincides with the decay rate dictated by 
Eqs. CJ-ljSll- This is the so-called on-shell G'^W approx- 
imation, which within RPA is usually called the on-shell 
G^W'^ approximation. There are various ways of go- 
ing beyond this approximation, by either normalizing the 



13 



0.9 r 



0.8- 



o o o o o O o 



T 

s 



o o o o o 



o o o 



0.7- 



0.6L 



e - / eV ^ 



FIG. 12: The real part of the renormahzation factor Z{ei) 
of Eq. Mm . versus Si — ef, as obtained from a first-principles 
evaluation of the G^W^ electron self-energy in Cu (open cir- 
cles) and Pd (crosses). 



electron energy, introducing short-range xc effects, or by 
looking at the so-caUed spectral function within a non- 
sclfconsistent or a selfconsistent scheme for the Green 
function. 



Electron-energy renormalization 

Near the energy shell, i.e., assuming that the deviation 
of the complex quasiparticle energy Ei from its noninter- 
acting counterpart Si is small, Ei can be expanded as 
follows 



E, w + AS(£,) + iE~e, 



9AI](w) 



This expansion yields 

E,^e^ + Z{e^)AJ:{e^), 



(39) 



(40) 



where Z is the complex quasiparticle weight or renormal- 
ization factor: 



1 - 



9AE(w) 



(41) 



RPA calculations of the renormalization factor Z{ei) 
of a FEG were performed by Hedin.j73| These calcula- 
tions show that in the case of a free gas of interacting 
electrons the renormalization factor near the Fermi level 
varies from unity in the high-density (r^ — > 0) limit to 
^ 0.6 at rg — 6, the imaginary part being very small. 
We have carried out first-principles RPA calculations of 
the renormalization factor of a variety of real solids. We 
have found that as in the case of a FEG at metallic den- 
sities (r^ = 2 — 6) this factor varies near the Fermi level 
from ^ 0.8 to ~ 0.6, as shown in Fig. 1121 where the real 



part of the renormalization factor Z{ei) of Cu and Pd is 
represented versus the noninteracting electron energy e^. 
Introducing Eq. into Eq. one finds 

r^i = -2 [ReZ(e,)ImAI](e,) -I- ImZ(£,)ReAS(e,)] . 

(42) 

Since the imaginary part of the renormalization factor 
Z{ei) is typically very small, Eq. 142() yields quasipar- 
ticle lifetimes that are larger than those obtained on- 
the- energy -shell by 30 — 40%. However, one must 
be cautious with the use of Eq. (|42|l when the self- 
energy is calculated in the G'^W approximation by re- 
placing the electron Green function by its noninteracting 
counterpart, which is built from the single-particle wave 
functions <?!>i(r) and energies In this approximation, it 
should be more appropriate to evaluate the self-energy 
on-shell by also replacing Ei by St in the calculation of 
the self-energy. This on-shell approximation yields a life- 
time broadening of the form of Eq. (|42|l with Zisi) = 1, 
which exactly coincides with the decay rate dictated by 
Eqs. CJ-©. 



GWT approximation 

Exchange and short-range correlation of the excited 
hot electron with the Fermi system, which are absent in 
the G'^W approximation, can be included in the frame- 
work of the GWT approximation. .74. ,J5] In this approx- 
imation, the electron self-energy is of the GW form, 
but with the actual time-ordered screened interaction 
W-^'^(r, r'; w) replaced by an effective screened interac- 
tion W'^'-'{y, r'; Lo) which for > is obtained as follows 



VF(r, r'; w) = w(r, r') -I- j dridv2 ['^(r, ri) 
+ fxc{-r,Vi]uj)]x{vi,V2]Uj)v{v2,v'), (43) 
the density-response function x(r,r';a;) now being 

X(r,r';w) = x*'(r,r';w) j dvi j dr2 x°(r,ri;w) 
X [t;(ri,r2) + /^c(r,ri;w)] x(r2,r';cj). (44) 

The kernel fxc{j^, ri; w) entering Eqs. H43|l and H44|l . which 
equals the second functional derivative of the xc en- 
ergy functional Ej;c[n{r)], accounts for the reduction in 
the electron-electron interaction due to the existence of 
short-range xc effects associated with the excited quasi- 
particle and with screening electrons, respectively. In 
the so-called time-dependent local-density approximation 
(TDLDA)'?^! or, equivalently, adiabatic local-density ap- 
proximation (ALDA), the exact kernel is replaced by 



ALDA 



(r,r';c.) 



[nebcin)] 



dn? 



5{r 



(45) 



n— n(r) 
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where Excin) is the xc energy per particle of a uniform 
electron gas of density n, and n(r) is the actual electron 
density at point r. 

Introduction of the GWT self-energy into Eq. H38|) 
yields on-the-energy-shell {Ei Si) an expression for 
the decay rate of the form of Eqs. Q"® with the 
screened interaction W{r,r';uj) of Eq. ^ replaced by 
that of Eq. (|^ . Existing first-principles calculations 
of hot-electron lifetimes in solids have all been carried 
out by neglecting short-range xc effects, i.e., by taking 
/2;c(r,r';w) = in Eos. and El First-principles cal- 
culations with full inclusion of short-range xc effects have 
been performed only very recentlv. (tt)] 

Spectral function 



4000 1 r 




8. - e„/eV -> 

I F 



In the framework of many-body theory, the propaga- 
tion and damping of an excited electron (quasiparticle) 
in the one-particle state (j>i{r) are dictated by the peaks 
in the spectral function Ai{uj), which is closely related 
to the imaginary part of the one-particle Green function 
G(r,r';c.):72j 

A{iu) = ^JdrJ dr'cb*{v)lmG{r,r'-iuWr'), (46) 
where 

- [E(r,r';tj)-y(r)5(r,r')]. (47) 

The noninteracting Green function Go(r,r';a;) is ob- 
tained from the eigenfunctions [cjji (r)] and eigenvalues [si] 
of an effective single-particle hamiltonian [H ~ — Vr/2 + 
Vir)]. 

The energetic position of the peak in the spectral func- 
tion defines the real part of the quasiparticle energy. The 
damping rate or lifetime broadening of the quasiparticle 
is given by the full width at half maximum (FWHM) 
of the peak: 

r,-' - A,. (48) 

Assuming that the peak in the spectral function Ai{LL!) 
has a symmetric Lorentzian form, this definition of the 
damping rate coincides with the decay rate dictated by 
Eq. jSHJ. 

We have carried out first- principles G°VF° calculations 
of the spectral function in Al. For an excited electron 
at the W point with energy Si — 1.02eV above the 
Fermi level, we have found the spectral function repre- 
sented in Fig. El This figure shows that the spectral 
function has indeed a single well-defined quasiparticle 
peak at the quasiparticle energy Sqp = 0.92 eV above 
the Fermi level with the FWHM = 24meV corre- 
sponding to the quasiparticle lifetime r = 27 fs. Similar 



FIG. 13: Spectral function of an excited electron at the W 
point in Al {si — 1.02 eV), versus Si — ef, as obtained from 
Eq. II46II in the G°W° approximation. 

calculations were performed by Schone et aL|l9l| for vari- 
ous one-particle Bloch excited states in the simple metal 
Al and the noble metals Cu, Ag, and Au. These authors 
found hot-electron lifetimes that are in reasonable agree- 
ment with the G°W° calculations obtained by Campillo 
et aZ.QIlSllll from Eq. 

Full self-consistent GW self-energy calculations, where 
the interacting Green function G(r, r';cj) of Eq. H47|l is 
used self-consistently to evaluate both the screened in- 
teraction and the GW self-energy, have been carried out 
recently for the FEgItII and simple semiconductors [tqI 
Is^l. Although self-consistency is nowadays known to 
yield very accurate total energies, existing calcula- 
tions tend to indicate that non-self-consistent calcula- 
tions should be preferred for the study of quasiparti- 
cle dynamics, [t^ [t^ This is due to the fact that ver- 
tex corrections not included in the self-consistent GW 
self-energy might cancel out the effect of self-consistency, 
thereby full self-consistent self-energy calculations that 
go beyond the GW approximation yielding results that 
might be close to G'^W'^ calculations. 



SUMMARY AND CONCLUSIONS 

We have presented a survey of current theoretical in- 
vestigations of the ultrafast electron dynamics in metals. 

First of all, a theoretical description of the finite inelas- 
tic lifetime of excited hot electrons in a many-electron 
system has been outlined, in the framework of time- 
dependent perturbation theory. Then, the main factors 
that determine the decay of excited states have been dis- 
cussed, and the existing first-principles theoretical inves- 
tigations of the lifetime of hot electrons in the bulk of 
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a variety of metals have been reviewed and compared 
to the available experimental data. The decay rates 
obtained within first-order time-dependent perturbation 
theory have been shown to coincide with the on-shell 
G'^W'^ approximation of many-body theory and to pro- 
vide a suitable framework to explain most of the available 
experimental data for simple and noble metals. 

Finally, various ways of going beyond the G'^W^ ap- 
proximation have been discussed, by normalizing the 
electron energy, introducing short-range xc effects, or 
looking at either the G°W^° or the self-consistent GW 
spectral function. First-principles GWT calculations 
with full inclusion of short-range xc effects have been 
carried out only very recently. [22| 

Alternatively, nonperturbative treatments of the inter- 
action of external excited electrons with a Fermi system 
have been developed. [^H^l These treatments are based 
on phase-shift calculations from kinetic theory[s3 and 
a modification of the Schwinger variational principle of 
scattering theory, 's^l both implemented for electrons in 
a FEG. The role of spin fluctuations in the screening and 
the scattering of excited electrons in a FEG has also been 
discussed'SS^ in the framework of kinetic theory, with the 
use of simple physically motivated models. 

The experiments considered here all involve very low 
densities of excited electrons, and the mutual interaction 
of excited electrons have been completely neglected in 
the theoretical calculations. In a situation in which the 
density of excited electrons is high, a different theoretical 
approach would be needed. Calculations al ong these lines 
have been carried out by Knorren et al. . |84L Isjj j by solving 
the Boltzmann equation for carriers in the conduction 
band. 
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